Practical techniques of analytic perturbation theory of QCD 

B.A. Magradze a 
Abstract 

The Lambert-W explicit solutions to the QCD renormalization group (RG) 
equation are considered up to fourth order in the MS scheme. We compare, 
systematically, these solutions with the conventional asymptotical (iterative) 
approximations and with the exact numerical solutions to the RG equation 
in the domain with three quark flavours. Applications of these solutions in 
analytic perturbation theory (APT) are discussed. Using these (Lambert- 
W, asymptotical and exact numerical) solutions we reconstruct the expansion 
functions for the non-power APT series in the space- and time-like regions. 
These expansion functions are examined in the infrared region. It is shown 
' that the Lambert-W solutions provide the excellent accuracy. 

o ' 
o 

q ■ 1 Introduction 

<u : 

Recently, in works [I] |2| the exact explicit 2-loop solution to the QCD renormalization group (RG) 
CN . equation was considered. This solution is given by 



. / denotes the number of active quark flavours at the energy scale Q, 3^ is the kth order coefficient 
of the QCD /^-function, b\ = /3i/(3q, A = A-jg^ is the conventional QCD parameter and W(C) is the 
Lambert W function. The branches of W are denoted Wk(C)i k = 0, ±1, . . . . An exhaustive review 
of the Lambert-W function may be found in Ref. j3J. In our notations the RG equation has the form 

£ : Q 2 das E J) = P f MQ 2 , /)) = -£ /fe +2 (Q 2 , /), (2) 

r- 1 

the normalization condition is a s (fi 2 ,f) = g 2 /(47r), where n is the renormalization point and g is 
the gauge coupling of QCD. In the class of schemes where the beta-function is mass independent 3 Q 
and Bi are universal 



and Q{ 



the results for 3^ (I] and /3g [S] in the modified MS (MS) scheme are 



?/ _ I /2857 5033/ 325/ 2 
(4vr) 3 18 - + 54 



% = 71Z*l— 7^ + ^H> W 



of 1 / 149753 neoA , / 1078361 6508 \. / 50065 6472. \, , 1093 , 3 

(5) 

Here £ is the Riemann zeta-function (£3 = 1.202056903...). 

Using (^Q) the analytical continuation of the coupling in the complex Q 2 plane has been performed 
[DlZHn] ( see a ^ so papers 0E1)- Afterwards, in paper 0|, a higher order running coupling (in arbitrary 
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MS-like renormalization scheme) was expanded in powers of the scheme independent explicit two- 
loop coupling (J2) c 

oo 

«^ 2 ) = £c«ap(Q 2 X (6) 

n=l 

then, a new method for reducing the scheme ambiguity for the QCD observables was proposed. A 
similar expansion (motivated differently for an observable depending on a single scale) has been 
suggested in [TU] . 

The exact solution (^Q) and the approximation © proved to be useful in the dispersive (analytic) 
approach to perturbative QCD. The dispersive approach ^2 EE21 EH EEH El EESj has been devised to 
extend the appropriately modified perturbation theory calculations towards the low-energy frontier. 
In work ^S] a particular version of the analytic approach, the analytic perturbation theory (APT), 
was suggested. According to this approach the standard perturbative series in powers of the running 
coupling a s (Q 2 ) are replaced by a novel, non-power series. The individual terms of these series 
(determined by the functions A n (Q 2 )) provide the rigor momentum-plane analyticity properties for 
the expanded quantity itself. APT has several (phenomenological) advantages over the conventional 
perturbation theory even for moderate energies in the five-flavour region |17 | I18[ IT9]. A striking 
feature of the non-power expansions in APT is that they exhibit improved convergence properties as 
compared with the standard perturbation theory |15[ I18 ( [T9]. Thus, the large ^-contributions to the 
expansion coefficients of time-like observables, inherent in standard perturbation theory |22U23| . are 
automatically summed in the individual terms of the non-power series. For various phenomenological 
and theoretical consequences obtained within APT we refer to works jlHl I2()j. 

In works |S] and [23J the Euclidean and Minkowskian sets of functions d {A n (Q 2 )} and {2l n (s)} 
have been explicitly constructed. The Lambert-W solutions (^Q) and © were used. The results, up 
to third order, were given in a compact form. Thus, the integrals for the time-like set of functions 
(see below formula (|18jl ) were analytically performed. For both sets of the functions an elegant 
equations have been derived; in the 2-loop case [S] and afterwards, also, in all orders of perturbation 
theory [23]. These equations generalize the one loop equations obtained previously in works |1£ 



(k) 

Thus functions A n satisfy the equations 

d ~^M 1 = -»£ P f N Ai% +1 (Q% n = l,2.... (7) 

V N=0 

Global variants of these functions with the quark threshold effects have been also reconstructed 

In this paper we estimate accurateness of the explicit Lambert-W solutions © and the con- 
ventional asymptotical approximations (see below formula <|13jl) at low scales. We keep in mind 
applications of these solutions in conventional and analytic perturbation theories. As a standard 
for comparison the "exact" numerically determined running coupling is used. In sect. 2, we present 
relevant formulas for the running coupling up to 4-loop order and give numerical results obtained 
from the approximate and exact solutions. In sect. 3, we present theoretical formulas for the space- 
like and timelike "analytic" images of powers of the running coupling (the functions A n (Q 2 ) and 
2l n (s)). To reconstruct these functions, the Lambert-W as well as the asymptotical solutions are 
applied. In the infrared region we give numerical results obtained from these approximations to the 
images, together with the "exact" results obtained by using a numerical method. Sect. 4 contains 
conclusions. 



c here and hereof we omit the argument /. 

d i.e. the functions of the Euclidean and Minkowskian arguments. 
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2 The Lambert-W solutions vs the asymptotical solutions 

In practice, it is sufficient to retain the first few terms in the expansion (jHJ). Then the approximation 
to the coupling reads 

n=N tr 

«£(Q 2 ) = £ ^^Hq 2 ) (8) 

n=l 

where k = 3,4... and = 1. The coefficient is arbitrary; this reflects the arbitrariness in the 
definition of the A parameter. The conventional MS definition of A [25] corresponds to Cg = 0. The 

(k) 

coefficients ch , for n > 2, may be readily determined. In the 4-loop order, the first nine coefficients 



are 



,(4) _ 02 (4) _ 03 (4) _ 5 f3% X 0: 



2 ' 



(9) 



3 /V 4 2/3 ' 5 3 01 601 

(4) 1 0102 2 1 0301 2 ^0203 

12 /3 3 + 12 01 +I 01 ' 

(4 ) = }_0\g_ 16 /g _ 4ftftft 11 /g _ 1 /33/g? 
C? 20 /? 4 5 /3 3 5 01 20 /3 2 20 /3 4 ' 

(4) = 19 0301 _ 23 0101 9 020301 _ 4V_0i0l _ 1 0\0l ]_030f 
Cg 3 0l 60 /5 4 20 /3 4 120 0l 30 30 0° 

(4) _ 946 /?!/3 2 2 /3 3 41 010203 134, 0201 2069 /? 2 4 103/gffj 149/??/% 1 g/g 1 ffift 
° 9 315 /3 4 140 fl§ + 35 01 315 /3 4 420 0° 504 /3 4 42 /3g 42 0* ' 

The coefficients have been calculated with symbolic manipulation system Maple V (release 5). It is 
convenient to apply the system Maple V to carry out numerical calculations, when using analytical 
expressions for A n (Q 2 ) and 2l ra (s) determined in terms of the Lambert W- function (see below). 
Maple has an arbitrary precision implementation of all branches of the Lambert-W function. On 
the other hand, for numerical problems professional FORTRAN compiler is preferable. This allows 
one to reduce the necessary computer time. However, the Lambert-W function is not available 
in standard FORTRAN libraries. To avoid this obstacle, one may numerically integrate the RG 
equation in the complex Q 2 -plane, instead of using the Lambert-W solutions. To the kth order the 
implicit solution to the RG equation is given by 



ln (t^V) = - - fo i ln ( 6 i + -)+ r ( =^y J dx > ( 10 ) 

where we have introduced the quantities a s ( ® 2 ) = 0oa s (Q 2 ) and <yk \a s ) = 0o0^ k \if-), and 

bi = In the 3- and 4-loop cases the integral on the right (1101) can be done. In the 3- loop case we 

Po 

find 



In ( — ^ = ln ( — 7T H — - + 62^ H 2 1 ( arctan ( — 

UW7 a s 2 W s a s 2 J ^A© V V 



+ 2b 2 a s \ ( h 

arctan 



(11) 
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where b k = l3 k /Po +1 and A( 3 ) = 4b 2 - b\. Formula (JTTJ) is valid only for A( 3 ) > (in the MS scheme, 
for / = 3, A® « 2.91). In the 4-loop case the result is more complicated 

^ = ~ ~ h In (h + 1) + , J®{a.) 

a s V a sJ hO- + b 1 xi) 



62 + 63X1 ({a s -xi) 2 n\ b 2 - h/h , / K + l/bi) 2 bjn 



j(4) (0s) = ^T" 3 : 1 In i"'-,V - 17 d "\ In 



62 + 632:1, LO x 62-63/61/ 2\\ 1 / /2a s + m\ / m \ 

— — — r— (m + 2xi) — — — - jz , m — — ; arctan — . — arctan . 



(12) 



where we denote 

m = u + v + —r, n = -m 2 + -(u — i>) 2 , A^ = 3(u — t>) 2 , xi = u + v — —r, 

U{x) = x 2 + mx + n, n= (\/^-g/2) 1 / 3 , u = -(V* + g/2) 1/3 , 

\3 , / /0 \2 _ 1 2 _ 2 3 , 1 6 2 61 t /?fc 



<5 = (p/3) + (g/2) , p = s--r, q = —r--rs + —, r = —,s = —,b k 



3' ' ^"27' 3'° ' 63' "63' &3' ""/jff +1 " 

Formula l|12|) has been obtained under the assumptions <5 > and \/<5 > q/2 (in the MS scheme, for 
f = 3, y/5 - q/2 m 0.003). 

In work [27], the 3-loop Eq. Qll[) has been numerically solved in the complex Q 2 plane. In this way 
the spectral densities (discontinuities) were numerically calculated. Evidently, from the technical 
point of view, this approach is more complicated than the using of the Lambert-W solutions. Strictly 
speaking, mathematical conditions for using of the numerical method should be investigated, i.e. the 
analytical structure of the coupling must be established. In a mathematically rigor manner this task 
has been performed only in the 2-loop case ^[2JH)|. Nevertheless, in this article, we shall use this 
method as well, up to fourth order in APT. 

In QCD, practitioners use asymptotical (iterative) solution to the RG equation. It is presented 
in the form of the asymtotical series |2H] 

, n 2\ _ 1 01 lnL 

a.,„W ) - ^ ^ L2 

+ ^(| (to2i - taL - I) + l) (13 » 



where L = lnQ 2 /A^ . The first line of Eq. (|13|) includes the 1- and 2-loop contributions to the 
coupling, the second line is the 3-loop correction and the third line is 4-loop correction. In the 3-loop 
case, coupling Q13JI has been frequently used in APT for calculation several spacelike and timelike 
observables |27M29| . On the other hand, for small momentum transfer the approximant (|13|) differs, 
significantly, from the exact numerical solution to the RG equation (see Table 1). In fact, formula 
(|13jl is valid for L ^> 1. Therefore, application of (|13j) . within infrared region, requires a justification. 

Throughout in this paper we consider QCD with / = 3 active flavours. For the QCD parameter 
we take the reference value A^^ = 0.400 GeV uniformly for all approximants considered. One could 
use a more sophisticated method to compare the approximants, but this could have caused small 
changes of the final results. 

In Table 1, we compare the Lambert-W solution (JHJ) to third order (with Nt r = 5) e and the 
asymptotical solution (|13|) (to second and third orders) with the 3-loop "exact" numerical cou- 
pling a n um, the numerical solution of Eq. Diff.(%,) stands for the percentage deviations: 

e The 3-loop order coefficients c|? can be read from @ setting fa = 0. 
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Figure 1: Plot of the functions oinlm{Q 2 ) and aif(Q 2 ,N tr ) ( N tr = 4,6,9). The solid line refers to 
the "exact" numerical coupling anlm{Q 2 )- 

Diff.(%,ts.3)= (1 - OL { S/a { nl m ) * 100, Diff.(%,as.k)= (1 - a&V«num) * 100, k = 2,3. The best 

(3) 

accuracy is achieved with the Lambert-W approximant a\ s . 

In Table 2, we summarize the results obtained from the 3- and 4-loop couplings. a>num denotes 
the "exact" numerical solution of Eq. (jTU)) in the 4-loop case (see Eq. ((T2*|0 . a[f is the 4-loop 
Lambert-W solution (JHJ) for N tr = 9 and etas (where k=3,4) refers to the asymptotical coupling (fT3)) 
at the 3- and 4-loop orders. We see again, that the best accuracy achieved with ocf s \ Curiously, in 

(3) 

this region, the 3-loop asymptotical coupling a as gives a better approximation to the exact 4-loop 
coupling than a^J ■ 

In Table 3 we give the percentage deviations of the Lambert-W approximations @ to fourth 
order (for N r =4,6,9) from a4t)m- We see, that the approximation (JSJ) is gradually improved when 
Nt r increases. The convergence of the series ((BJ worsens only close to the Landau singularity (in the 
4-loop case the Landau singularity occurs at ~ 1.773Aig^ ~ 0.709 GeV for / = 3 and A-^r = 
0.400 GeV). In Figure 1 we plot the "exact" numerical coupling otnlm{Q 2 ) and the approximants 
aig (Q 2 , Nt r ) (for Ntr = 4, 6, 9) as a function of the momentum variable Q = \J Q 2 . 



3 The Euclidean and Minkowskian expansion functions 

For the reader's convenience here we briefly summarize main aspects of APT. Let D(Q 2 ) be Adler 
D-function related to some timelike process. It has the asymptotic expansion in powers of a running 
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Table 1: Percentage deviations of the approximants a ts (N tr = 5), a a J and a as from the exact 
(numerical) 3- loop coupling anlm- 



\fQ i GeV 


(3) 

C-t Yi 11 'til 


Diff.(%,ts.3) 


Diff.(%,as.3) 


Diff.(%,as.2) 


.80 


0.76491 


1.87 


-15.49 


-7.16 


.90 


0.63323 


0.81 


-9.25 


-3.93 


1.00 


0.55414 


0.43 


-6.08 


-1.59 


1.10 


0.50028 


0.26 


-4.33 


0.01 


1.20 


0.46075 


0.17 


-3.29 


1.12 


1.30 


0.43025 


0.12 


-2.62 


1.90 


1.40 


0.40587 


0.09 


-2.17 


2.47 


1.50 


0.38583 


0.06 


-1.86 


2.88 


1.60 


0.36901 


0.05 


-1.63 


3.19 


1.80 


0.34220 


0.03 


-1.33 


3.60 


2.00 


0.32165 


0.02 


-1.14 


3.84 


2.2 


0.30527 


0.02 


-1.01 


3.98 


2.6 


0.28059 


0.01 


-.85 


4.11 



coupling 

oo 

D P t(Q 2 ) =D (1 + J2 d n a n s {Q 2 )), (14) 



n=l 



here Dq is a process dependent constant. In APT, the same quantity should be presented in the 
form of a nonpower asymptotic expansion as 

oo 

D an (Q 2 ) = D (l + Y,d n A n (Q 2 )) (15) 

n=l 

where A n is the "analyticized nth power" of the coupling in the spacelike region. This is determined 
through the spectral representation 

A n {Q 2 ) = - ^ -^^da = - r t J > 2 , Pn(t)dt, (16) 
Wo cr + Q n J-oo (e + Q7 A ) 

where the spectral function is defined as 

p n (a) = 3?K(-<7-*0)}, 

A on the right of (fTo]) denotes the QCD scale parameter and p n (t) = p n (cr) with t = ln(^). Let R(s) 
be the physical quantity reconstructed through D(Q 2 ) in the timelike domain (notable examples are 
R e + e -( s ) = c(e + e~ — > hadrons) /a(e + e~ — > p. + u~) and R T (M T ) = F(r — ► hadrons) /T(t — > v T v s e). 
Then, in APT it has the representation [22] 



R(s) = R (l + r(s)) where r(s) = £^ =1 d n % n (s). (17) 
The timelike set of functions {2l n (s)} is defined by the formula fZT\ 

1 roc 1 

Onto = - / —Pn{cr). (18) 
7T J s a 

We see, that the spectral functions corresponding to powers of the coupling are of special importance. 
For < / < 6, the spectral function corresponding to the nth power of the coupling (^Q) reads [H1IH1 

pg'W = (|)"a(- 1 + ^ Nt) J ; z(t) = ^(-l- t , bl+ >« Wbl -l)), (19) 
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Table 2: The approximants a ts (N tr = 9), a a J and a as vs the 4- loop "exact" numerical coupling 
ttrmm- Percentage deviations of these approximants from oenlm are given. 



GeV 


(4) 


(4) 
a ts 


(4) 

a as 


(3) 
aas 


Diff.(%,ts) 


Diff.(%,as4) 


Diff.(%,as3) 


0.80 


0.91262 


0.86804 


1.04504 


0.88340 


4.89 


-14.5 


3.20 


0.90 


0.69243 


0.68439 


0.76687 


0.69179 


1.16 


-10.7 


0.09 


1.00 


0.58701 


0.58462 


0.63081 


0.58784 


0.41 


-7.46 


-0.14 


1.10 


0.52156 


0.52063 


0.54982 


0.52195 


0.18 


-5.42 


-0.07 


1.20 


0.47585 


0.47542 


0.49548 


0.47589 


0.09 


-4.12 


-0.01 


1.30 


0.44164 


0.44142 


0.45606 


0.44153 


0.05 


-3.26 


0.03 


1.40 


0.41484 


0.41471 


0.42589 


0.41469 


0.03 


-2.66 


0.04 


1.50 


0.39313 


0.39305 


0.40190 


0.39301 


0.02 


-2.23 


0.03 


1.60 


0.37510 


0.37505 


0.38224 


0.37503 


0.01 


-1.90 


0.02 


1.70 


0.35984 


0.35980 


0.36577 


0.35982 


0.01 


-1.65 


0.00 


1.80 


0.34670 


0.34667 


0.35173 


0.34674 


0.01 


-1.45 


-0.01 


2.00 


0.32515 


0.32514 


0.32892 


0.32530 


0.00 


-1.16 


-0.05 


2.20 


0.30811 


0.30810 


0.31106 


0.30835 


0.00 


-0.96 


-0.08 


2.40 


0.29421 


0.29421 


0.29660 


0.29452 


0.00 


-0.81 


-0.10 


2.60 


0.28261 


0.28261 


0.28459 


0.28297 


0.00 


-0.70 


-0.13 



one may rewrite the spectral function (|19|) in the equivalent form, in terms of the branch W—i, using 
the symmetry property Wk(z) = W_k{z) A- Thus, we find 

^--(t)M-ITOT)"- 

where z\(t) = z{t). The spectral function corresponding to the coupling (jHJ then is 

Ntr 

p ( fs\(t) = E c( n ) P { n ) (^ (20) 
n=l 

this gives the following representations for the first expansion functions 

n=N tr n=N tr 

-4 fc) i(Q 2 ) = E ^AtUQ 2 ) ^d <\( S ) = £ cW< n (s), 

n=l n=l 

analogical formulae hold for the functions with a higher index. 

In practice integral (|16|) should be regulated. In Ref 8] the following formula was derived 

a to 2 ) - a (n 2 m + 21 ( A 2 ,^ + / 0{Q 2 e- R /R l+n ) if Q 2 > 1 

An{Q ) - An{Q ,K) +^n(A e ) + < 0(Q~ 2 e~ R /R 1+n ) if Q 2 < 1 ^ ' 

where A n (Q 2 ,R) denotes the integral (|16|) taken over the finite interval —R\ < t < R, functions 
2t n , n = 1,2..., are defined in (|18|) and Q 2 = Q 2 /A 2 . For sufficiently large values of R, when 
Q 2 exp(— R)/R 1+n <C 1, the contributions of order e~ R can be omitted. So that the extra term 
2t n (A 2 e R ) compensates main error emerging because of truncation of the integral on the upper 
bound. Formula (|21|) enables us to achieve a good numerical precision even for moderate values of 
the cutoff R. In the region t — ► — oo the integrand decreases rapidly, so that we can take Ri < R. 
To achieve a good numerical precision it suffices to limit the integral to the region — 20 < t < 1000. 
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Table 3: Relative errors of the 4-loop approximant a ts (Q 2 , Nf S ), with Nt r = 4, 6, 9. The approximant 
is compared with the exact numerical 4-loop coupling oinlm- 



To 1 GeV 


(4) 


Diff.(%,ts) 

N tr = 4 


Diff.(%,ts) 

Ntr = 6 


Diff.(%,ts) 

N tr = 9 


.80 


.91262 


15.7 


9.4 


4.9 


.90 


.69243 


7.4 


3.3 


1.2 


1.00 


.58701 


4.3 


1.6 


0.4 


1.10 


.52156 


2.9 


0.9 


0.2 


1.20 


.47585 


2.1 


0.5 


0.1 


1.30 


.44164 


1.6 


0.4 


0.1 


1.40 


.41484 


1.2 


0.3 


0.0 


1.50 


.39313 


1.0 


0.2 


0.0 


1.60 


.37510 


0.8 


0.1 


0.0 


1.70 


.35984 


0.7 


0.1 


0.0 


1.80 


.34670 


0.6 


0.1 


0.0 


1.90 


.33525 


0.5 


0.1 


0.0 


2.00 


.32515 


0.5 


0.1 


0.0 



To calculate the Minkowskian functions (|18|) one may use a method of integrating expressions 
containing inverse functions Using this technique, in work [H] (see also papers |24l I30j ^ the 
following expressions have been obtained 

2lf) (s) = + J_ Im f 1 



Pi ttA V« (2) (-s), 
«?>M = -^/ m ( ta (l + |a<*><-.))), (22, 



where 



s + l 0) = -^.(l + W„ 1 {z s ))- 1 with z s = 6^ 1 (sA- 2 )- 1 / fel exp(-l-zvr(6^ 1 -l)). 
Pi 

The Minkowskian functions with a higher index are easily determined by using the recursion formula 

^xM( s ) = -n(m%(s) + (3i%%( S )) n = l,2..., (23) 
obtained in paper |Sj. 

So far as the asymptotical solution (|13|) is concerned one may use the following formula 

L~ n ln m L = (-l) m -^-L" n , (24) 
where L = L(Q 2 ) = ln(Q 2 /A 2 ), m and n are positive numbers. We can write, for a > 0, 



L(-a - iO) = R a exp (z<3?(cr)), where R a = \/ln 2 a + 7r 2 with o- = — , (25) 

and 

- arcsin {n / R a ) if a > 1, 



$(<r) 



-7r + arcsin (tt/R^) if < a < 1. 
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Table 4: The second order Euclidean expansion functions A n (Q 2 ) and A a s,n{Q 2 ), for n = 1,2, as 
a function of the momentum transfer yj Q 2 . A^ = 3 = 0.400 GeV. 



TO 1 GeV 


A? 


A [2) 

as 1 


Diff.(%,l) 


A {2) 


A (2) 


Diff.(%,2) 


0.4 


0.507853 


0.476767 


6.1 


0.118913 


0.113524 


4.5 


.60 


0.438444 


0.413090 


5.8 


0.107177 


0.100299 


6.4 


.80 


0.393408 


0.372388 


5.3 


0.09703 


0.090135 


7.1 


1.00 


0.361380 


0.343505 


5.0 


0.088705 


0.082211 


7.3 


1.20 


0.337219 


0.321662 


4.6 


0.081875 


0.075872 


7.3 


1.40 


0.318218 


0.304414 


4.3 


0.076208 


0.070681 


7.3 


1.60 


0.302807 


0.290363 


4.1 


0.071443 


0.066345 


7.1 


1.80 


0.290003 


0.278638 


3.9 


0.067385 


0.062665 


7.0 


2.00 


0.279159 


0.268669 


3.8 


0.063886 


0.059498 


6.9 


2.20 


0.269831 


0.260062 


3.6 


0.060839 


0.05674 


6.7 


2.40 


0.261702 


0.252536 


3.5 


0.05816 


0.054315 


6.6 


2.60 


0.254538 


0.245885 


3.4 


0.055784 


0.052163 


6.5 



Table 5: The second order Minkowskian expansion functions 2tn (s) and 2lL,n(s) (for n = 1,2) as a 
function of the energy scale \fs. Aj = 3 = 0.400 GeV . 



V* GeV 


0™ 


0,(2) 


Diff.(%,l) 


a?' 




Diff.(%,2) 


0.4 


0.495247 


0.458395 


7.44 


0.133488 


0.120209 


9.9 


.60 


0.417309 


0.393376 


5.73 


0.114407 


0.102211 


10.7 


.80 


0.369053 


0.352141 


4.58 


0.098808 


0.089398 


9.5 


1.00 


0.336208 


0.323077 


3.91 


0.087040 


0.079716 


8.4 


1.20 


0.312294 


0.301367 


3.50 


0.078109 


0.072190 


7.6 


1.40 


0.294001 


0.284473 


3.24 


0.071173 


0.066212 


7.0 


1.60 


0.279480 


0.270906 


3.07 


0.065653 


061368 


6.5 


1.80 


0.267616 


0.259734 


2.95 


0.061160 


0.057372 


6.2 


2.00 


0.257700 


0.250344 


2.85 


0.057432 


0.054022 


5.9 


2.20 


0.249258 


0.242318 


2.78 


0.054286 


0.051172 


5.7 


2.40 


0.241962 


0.235360 


2.73 


0.051594 


0.048718 


5.6 


2.60 


0.235575 


0.229257 


2.68 


0.049261 


0.046581 


5.4 



Evidently, the function L~ n ln m L for m, n > satisfies the Kallen-Lehmann representation. Let 
p(a, m, n) be the corresponding spectral density. Combining (|24j) and (|25|) we readily derive the 
following formula 

p(a,m,n) = (-l) m+1 (i?-" s in(n$(a))) . (26) 

Formula (|26|) allows one to derive explicit expressions for the spectral functions corresponding to 
the asymptotical coupling (fT3*|) (and to its positive powers) in various orders of perturbation theory. 
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Thus, we have 



(i) / \ _ 77 

Pas,l\ a ) — a /, 2 - , 9^ 

.(2) _ „(i) ^ A 



7T fli 1 (27) 

fl n 2 - ~ 2 , " % n 2 , 2 (vrlnaln(ln 2 a + ^ 2 ) + (In 2 a - vr 2 )$(a)), 

/?o (In cr + vr 2 ) /3g (lnV + vr 2 ) 2 

p£iW = pSiW + ^ (||(pfo2,3) - p((7, 1,3)) + (| " |) PM,3)) , etc. 

To reconstruct the timelike analytic images of the powers of the coupling (|13|) we define the auxiliary 
quantities 

1 Z" 00 

2l(s,m,re) = — / p(a,m,n)da. (28) 

In the case s > A 2 , one may change the integration variable in (|28|) according it = arcsin(7r/iZ cr ). 
We thus find the representation 

a / /-arcsin (tt/R s ) \ 

2l(s,m,n) = — ( *~ n J sm n ~ 2 usm{nu)du\ , (29) 

here use has been made of (|26|). For < s < A 2 , the analogical formula can be written. In the 
most cases of practical interest the integrals (|29|) can be performed in terms of elementary functions. 
Using Q29|) . for example, in the 2-loop case, if s = s/A 2 > 1, we find 

<M = } o (0J> - Cretan (!=!)) + |i, (-1 - In*, ^arcsin (£)) , (30) 



where -R s = \/ln 2 s + tt 2 and s = s/A 2 

In Table 4 we present numerical results obtained with the Euclidean expansion functions in the 
2-loop case. An\Q 2 ) stands for the Euclidean analytic image of the nth power of the exact explicit 
2-loop coupling l(T|). The Euclidean analytic images of the powers of the 2-loop asymptotical coupling 
(the first line in Eq. (|T3|) ) are denoted as A^as,n(Q 2 ), n = 1,2 . . .. We examine the first two Euclidean 
functions (n = 1, 2) in the region with three flavours 0.4 GeV < \/Q 2 <2.6 GeV. We see, that the 
errors in the approximations A^ n (Q 2 ) for n = 1, 2 are large. The relative errors are determined as 

Diff.(%,n)= (1 - 4?n(Q 2 )M 2) (Q 2 )) * 100. 

In Table 5 we give the 2-loop order results obtained with the corresponding Minkowskian ex- 

(2) (2) 

pansion functions 2ln (s) and 2las, n (s) fo r n = 1,2 in the region 0.4 GeV < < 2.6 GeV. It is 

(2) (2) 

seen from the table that the errors (Diff.(%,n)= (1 - W as >, n (s)/K >(s)) * 100) in the asymptotical 

(2) 

approximants 2las,n(s) are too large. 

In Tables 6 and 7 we estimate the accuracy of the 3- loop order approximants A^ n (Q 2 ) and 
2Q S n (s) that are constructed from the truncated series (jHJ) where Nt r = 5. In the domain with three 

flavours, we compare these functions with the "exact" functions Anlm,n{Q 2 ) and %lnlm,n{s). The 
"exact" functions are reconstructed numerically by solving the transcendental equation (|lll) in the 
complex-Q 2 plane. We see, that both (Euclidean and Minkowskian) Lambert-W approximants give 
the excellent accuracy in the considered region. Beyond this region, these approximations are even 
more accurate. 

In Table 8, we summarize the numerical results obtained with the 3-loop Euclidean functions 
AfniQ 2 ) and A { aln{Q 2 ) for n = 1, 2. We take again N tr = 5. The functions A^ n (Q 2 ) correspond to 
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the 3-loop asymptotical coupling. We see that the theoretical errors (Diff.(%,n)= (1 — A a s,n/ A\ s n ) * 

100) in A [ aln{Q 2 ) for n = 1,2 are significantly smaller than those in the 2-loop case. The 3-loop 
results for the Minkowskian expansion functions are given in Table 9. We see, that the errors in 

(3) 

the approximants %las\n{s) are approximately of the same size as those in the corresponding 3-loop 
Euclidean expansion functions. 

In Tables 10 and 11, the results obtained in the 4- loop case for the Euclidean and Minkowskian 
expansion functions are given. A^ m 1 (Q 2 ) and 2l^ 4 2 m i( s ) stand for the "exact" numerical functions 
obtained by solving the transcendental equation (|10j) (to fourth order) in the complex Q 2 -plane 
(see Eq. ©). The functions J$\(Q 2 ) and %!j£\(s) are reconstructed from the truncated series (181) 
where N tr = 6. The approximants A^ s \(Q 2 ) and 2ljj^i(s) correspond to the asymptotical coupling 
(fT3*|) . We reproduce practically exact results with the approximants A^\(Q 2 ) and 2l£j\(s). The 

asymptotical approximants A^Ji(Q 2 ) and 21^1 ( s ) ^ ea< ^ to the results, also, in fairly good agreement 
with the exact functions. 

4 Conclusion 

The explicit Lambert-W solutions and the conventional asymptotical approximations to the QCD RG 
equation are compared with one another in the MS scheme up to 4- loop order. These approximations 
have been carefully examined in the infrared region, where they are expected to be less accurate. As 
a standard for the comparison we have used the "exact" numerically calculated running coupling. 
Beyond second order, the Lambert-W approximations were represented as series in powers of the 
exact 2-loop coupling. We have shown, that the partial sums of these series, with the first few terms, 
give sufficiently good successive approximations to the "exact" numerical coupling. The required 
accuracy is achieved by choice of a sufficiently large value for the truncation order Nt r , even close 
to the Landau singularity (see Figure 1). This was confirmed to third, as well to fourth orders (see 
Tables 1 and 2) f . Note that, at the energy scale of the r-lepton mass M T = 1.77 GeV, for A = 0.400 
GeV, the error in ctfj (Q 2 ,N ts = 9) is about 0.01%, while the error in ahs (Q 2 ), at the same scale, 
is more sizable; Diff.(%,as)=1.6. 

Applications of these solutions in APT are discussed. We have presented convenient theoretical 
expressions for the Euclidean and Minkowskian analytic images of powers of the running coupling, 
which are determined in terms of the Lambert-W function (see sect. 3). Alternative expressions have 
been also derived using asymptotical solution (|13jl. In the 2-loop case, the asymptotical images 
(i.e. the images of the asymptotical coupling (|13|0 are found to lead to large theoretical errors. For 
example, the errors in the functions 21^ (m 2 ) and %^{M 2 ) (m p = 0.770 GeV and M T = 1.77 GeV 
being masses of the p mezon and T-lepton) are found to be 4.6% and 3% respectively. However, 
we have observed, to third and fourth orders, that the accuracy of the asymptotical Euclidean and 
Minkowskian approximants is significantly improved. Thus, the errors in for s = m 2 p and 

s = Ml are found to be 0.9% and 0.6% respectively (see Table 9). In the 4-loop case, these errors 
are even smaller; 0.4% and 0.5% respectively (Table 11). 

On the other hand, we have reproduced practically exact results with the Euclidean and Minkowskian 
images for the Lambert-W coupling (jHJ) to third and fourth orders (see Tables 8,9,10 and 11). An- 
other argument in favor to the Lambert-W solutions is that they enable one to derive closed analytic 
expressions for the time-like observables |241 130j . Of considerable importance in numerical and an- 
alytical calculations are the recursion formulas ® and (|23|) which enable one to write compact and 
clear Maple programs. 

f We shall give a mathematical investigation of the convergence properties of the series © elsewhere. 
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Table 6: The Euclidean approximant A^iiQ 2 ) reconstructed from the truncated series (JHJ) with 
Ntr = 5 vs the "exact" (numerical) 3-loop function A^l m ^{Q 2 ). A /=3 = 0.400 GeV. 



Ge y 


„4 (3) 

~^num,l 




Diff.(%) 


TO 1 GeV 


„4 (3) 




Diff.(%) 


0.4 


0.511938 


0.511933 


0.001 


1.8 


0.295068 


0.295069 


-0.000 


0.6 


0.443765 


0.443761 


0.001 


2.0 


0.284026 


0.284028 


-0.001 


0.8 


0.399121 


0.399119 


0.000 


2.2 


0.274510 


0.274511 


-0.000 


1.0 


0.367134 


0.367133 


0.000 


2.4 


0.266203 


0.266205 


-0.001 


1.5 


0.315504 


0.315505 


-0.000 


2.6 


0.258874 


0.258876 


-0.001 



Table 7: The Minkowskian approximant 2Lj s \(s) corresponding to the truncated series (JHJ where 
Nt r = 5 vs the "exact" (numerical) 3-loop function 9$lm,v A /=3 = °- 400 GeV ■ 



y/a GeV 


n,(3) 

num,l 


o,(3) 


Diff(%) 


V« GeV 


or(3) 

n«m,l 


o,(3) 


Diff(%) 


0.4 


0.501635 


0.501609 


0.005 


1.8 


0.273008 


.0.273011 


-0.001 


0.6 


0.425402 


0.425397 


0.001 


2.0 


0.262718 


0.262720 


-0.001 


0.8 


0.377104 


0.377110 


-0.002 


2.2 


0.253954 


0.253955 


-0.000 


1.0 


0.343725 


0.343733 


-0.002 


2.4 


0.246377 


0.246378 


-0.000 


1.5 


0.292420 


0.292425 


-0.002 


2.6 


0.239746 


0.239747 


-0.000 
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Table 8: The third order Euclidean approximants A ts n {Q 2 ) and A as \n(Q 2 ) (n=l,2) as a function 
of the momentum transfer ^Q 1 . Diff.(%,n)= (1 - A^, n /Af,n) * 10 °- 



GeV 


^ (3) 


A {3 \ 


Diff.(%,l) 


^ (3) 


.4 {3) 


Diff.(%,2) 


0.10 


754824 

v / . | '.J J. v 7 ^ -L 


754546 

KJ * 1 " 7 -L " X -L V / 


0.04 


0.118112 


115945 

V / • X X ".7 7 -L " 7 


1.83 


0.20 


635380 


638860 

VJ>V7UV7L7V7V_/ 


-0.55 


124962 


123223 

V / • x_ — ■ — — f X 


1.39 


0.30 


562731 


567327 


-0.82 


0.122311 


0.121454 


0.70 


0.40 


511933 


516651 


-0.92 


117273 


117054 


0.19 


0.50 


473770 


478275 


-0.95 


0.111773 


0.111952 


-0.16 


0.60 


443761 


447960 

V / • -L J- 1 ■ 7 V7 v 1 


-0.95 


106436 

V 7 . A. V / V / I * 7 V/ 


106854 

V / • X v 7 KJ \.J .J J- 


-0.39 


0.70 


419392 


423279 


-0.93 


101465 

V 7 . 1 V 7 X X V 7 '.J 


0.102025 


-0.55 


0.80 


399119 


402718 


-0.90 


096911 


097551 


-0.66 


0.90 


381931 


385274 


-0.88 


092763 

\J t\J ij i-l 1 Ky* } 


093446 


-0.74 


1.00 


0.367133 


0.370252 


-0.85 


0.088990 


0.089694 


-0.79 


1.20 


0.342862 


0.345611 


-0.80 


0.082421 


0.083127 


-0.86 


1.40 


0.323686 


0.326151 


-0.76 


0.076922 


0.077607 


-0.89 


1.60 


0.308076 


0.310316 


-0.73 


0.072264 


0.072920 


-0.91 


1.80 


0.295069 


0.297128 


-0.70 


0.068273 


0.068898 


-0.92 


2.00 


0.284028 


0.285938 


-0.67 


0.064815 


0.065410 


-0.92 


2.20 


0.274511 


0.276297 


-0.65 


0.061790 


0.062356 


-0.92 


2.40 


0.266205 


0.267885 


-0.63 


0.059121 


0.059661 


-0.91 


2.60 


0.258876 


0.260465 


-0.61 


0.056747 


.057262 


-0.91 



Table 9: The third order Minkowskian approximants W ts n (s) and 2lL,n(s) (n=l,2), in the domain 
with three flavours. Diff(%,n)= (1 - Slgn/St^J * 100. 



yfs GeV 


«(3) 
**t s 1 


n,(3) 


Diff.(%,l) 


«(3) 


n,(3) 

^as,2 


Diff.(%,2) 


0.10 


0.769627 


0.767795 


0.24 


0.126247 


0.122167 


3.23 


0.20 


0.640566 


0.648163 


-1.19 


0.140034 


0.136232 


2.71 


0.30 


0.559089 


0.568235 


-1.64 


0.138017 


0.137921 


0.07 


0.40 


0.501609 


0.508946 


-1.46 


0.130889 


0.132644 


-1.34 


0.50 


0.458689 


0.464136 


-1.19 


0.122465 


0.124524 


-1.68 


0.60 


0.425397 


0.429557 


-0.98 


0.114203 


0.116042 


-1.61 


0.70 


0.398821 


0.402184 


-0.84 


0.106619 


0.108168 


-1.45 


0.80 


0.377110 


0.379979 


-0.76 


0.099846 


0.101157 


-1.31 


0.90 


0.359032 


0.361579 


-0.71 


0.093863 


0.094997 


-1.21 


1.00 


0.343733 


0.346057 


-0.68 


0.088593 


0.089597 


-1.13 


1.20 


0.319204 


0.321235 


-0.64 


0.079832 


0.080667 


-1.05 


1.40 


0.300337 


0.302175 


-0.61 


0.072913 


0.073643 


-1.00 


1.60 


0.285311 


0.287003 


-0.59 


0.067343 


0.067999 


-0.97 


1.80 


0.273011 


0.274587 


-0.58 


0.062776 


0.063374 


-0.95 


2.00 


0.262720 


0.264199 


-0.56 


0.058965 


0.059517 


-0.94 


2.20 


0.253955 


0.255351 


-0.55 


0.055737 


0.056251 


-0.92 


2.40 


0.246378 


0.247702 


-0.54 


0.052967 


0.053448 


-0.91 


2.60 


0.239747 


0.241007 


-0.53 


0.050561 


0.051014 


-0.89 
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Table 10: The fourth order Euclidean approximants A\ S '^(Q 2 , N ts = 6) and A^J^Q 2 ) vs the "exact" 
numerical function A^ mjl {Q 2 ). A/= 3 = 0.400 GeV. 



y/Q 2 GeV 


A {4) i 


A {4) 




Diff.(%,ts) 


Diff.(%,as) 


0.1 


754123 


754128 


755096 


- 0.001 


-0.13 


0.2 


634500 


634501 


632852 


-0.000 


0.26 


0.3 


561 976 


561972 


559131 


0.001 


0.51 


0.4 


511361 


511355 


508229 


0.001 


0.61 


0.5 


473372 


473365 


470371 


0.001 


0.63 


0.6 


443510 


443504 


440810 


0.001 


0.61 


0.7 


419262 


419257 


416912 


0.001 


0.56 


0.8 


399087 


399083 


397085 


0.001 


0.50 


0.9 


381977 


381975 


380298 


0.000 


0.44 


1.0 


0.367243 


0.367242 


0.365854 


0.000 


0.38 


1.2 


0.343063 


0.343063 


0.342154 


-0.000 


0.27 


1.4 


0.323946 


0.323948 


0.323404 


-0.001 


0.17 


1.6 


0.308374 


0.308376 


0.308108 


-0.001 


0.09 


1.8 


0.295391 


0.295394 


0.295334 


-0.001 


0.02 


2.0 


0.284365 


0.284367 


0.284466 


-0.001 


-0.04 


2.2 


0.274857 


0.274859 


0.275080 


-0.001 


-0.08 


2.4 


0.266554 


0.266556 


0.266870 


-0.001 


-0.12 


2.6 


0.259226 


0.259228 


0.259613 


-0.001 


-0.15 



Table 11: The fourth order Minkowskian approximants 2lJ s l (s, N ts = 6) and % as \i(s) vs the "exact" 
numerical function 2l^ m ^s). Af =3 = 0.400 GeV. 



y/8 GeV 


ntim,l 


2l (4) 


%as,l 


Diff.(%,ts) 


Diff.(%,as) 


0.1 


0.768748 


0.768760 


0.772631 


-0.002 


-0.51 


0.2 


0.639012 


0.639034 


0.635957 


-0.003 


0.48 


0.3 


0.557764 


0.557745 


0.550459 


-0.003 


1.31 


0.4 


0.500760 


0.500717 


0.493308 


-0.009 


1.49 


0.5 


0.458283 


0.458247 


0.452271 


-0.008 


1.31 


0.6 


0.425335 


0.425315 


0.421015 


-0.005 


1.02 


0.7 


0.399005 


0.399001 


0.396164 


-0.001 


0.71 


0.8 


0.377464 


0.377470 


0.375788 


-0.002 


0.44 


0.9 


0.359498 


0.359509 


0.358695 


-0.003 


0.22 


1.0 


0.344270 


0.344284 


0.344107 


-0.004 


0.05 


1.2 


0.319808 


0.319822 


0.320436 


-0.004 


-0.20 


1.4 


0.300954 


0.300966 


0.301979 


-0.004 


-0.34 


1.6 


0.285916 


0.285925 


0.287128 


-0.003 


-0.42 


1.8 


0.273593 


0.273600 


0.274881 


-0.003 


-0.47 


2.0 


0.263275 


0.263281 


0.264579 


-0.002 


-0.50 


2.2 


0.254483 


0.254487 


0.255770 


-0.002 


-0.51 


2.4 


0.246879 


0.246882 


0.248133 


-0.001 


-0.51 


2.6 


0.240222 


0.240224 


0.241433 


-0.001 


-0.50 
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